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Abstract. Stein's method is used to study the trace of a random el- 
ement from a compact Lie group or symmetric space. Central limit 
theorems are proved using very little information: character values on 
a single element and the decomposition of the square of the trace into 
irreducible components. This is illustrated for Lie groups of classical 
type and Dyson's circular ensembles. The approach in this paper will 
be useful for the study of higher dimensional characters, where normal 
approximations need not hold. 



1. Introduction 

There is a large literature on the traces of random elements of compact 
Lie groups. One of the earliest results is due to Diaconis and Shahshahani 
[DS] . Using the method of moments, they show that if g is random from 
the Haar measure of the unitary group U(n, C), and Z = X + iY is a stan- 
dard complex normal with X and Y independent, mean and variance \ 
normal variables, then for j = 1,2, Tr{g 3 ) are independent and dis- 
tributed as y/]Z asymptotically as n — > oo. They give similar results for 
the orthogonal group 0(n,M) and the group of unitary symplectic matrices 
U Sp(2n,C). Concerning the error in the normal approximation, Diaconis 
conjectured that for fixed j, it decreases exponentially or even subexponen- 
tially in n. Stein [St2j uses "Stein's method" to show that Tr(g k ) on 0(n,M.) 
is asymptotically normal with error 0(n~ r ) for any fixed r. Johansson [J] 
proved Diaconis' conjecture for classical compact Lie groups using Toeplitz 
determinants and a very detailed analysis of characteristic functions. 

One direction in which the results of the previous paragraph have been 
extended is the study of linear statistics of eigenvalues: see [Jj, [DEj . [So] 
and the numerous references therein. There is also work by D'Aristotile, 
Diaconis, and Newman [DDN] on central limit theorems for linear functions 
such as Tr(Ag) where A is a fixed nx n real matrix and g is from the Haar 
measure of 0(n, R). In recent work, Meckes |Mc2] refined Stein's technique 
from |St2] to establish a sharp total variation distance error term (order 
n~ l ) for the [DDNj result. 
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A natural goal is to prove limit theorems (with error terms) for the dis- 
tribution of traces in other irreducible representations: i.e. x T (g)i where g 
is a random element of a compact Lie group and x T is the character of an 
irreducible representation r. This would have direct implications for Katz's 
work [Ka] on exponential sums; see Section 4.7 of |KLR| for details. We do 
not attain this goal, but make a useful contribution to it. 

More precisely, the current paper presents a formulation of Stein's method 
designed for the study of x r (<50- I n the case of normal approximation, we 
obtain 0(n~ l ) bounds for the error term using only two pieces of informa- 
tion: 

• The value of the "character ratios" J^^) wnere $ ma y be arbitrary 
but a is a single element of G (typically chosen to be close to the 
identity) 

• The decomposition of r 2 into irreducible representations 

In contrast, the method of moments approach requires knowing the mul- 
tiplicity of the trivial representation in r k for all k > 1 (which could be 
tricky to compute) and does not give an immediate bound on the error. 
Johansson's paper [J] gives sharper bounds when \ T is the trace of an el- 
ement from a classical compact Lie group, but required knowledge of high 
order moments and deep analytical tools which might not extend to arbi- 
trary representations r. Even the Stein's method approaches of Stein |St2] 
and Meckes [Me2| use information about the distribution of matrix entries; 
very little is known about this for arbitrary r, whereas the main ingredient 
for our approach (character theory) is well-developed. 

Let us explain our statement in the abstract that the methods of this pa- 
per will prove useful for approximation other than normal approximation. 
We use Stein's method of exchangeable pairs which involves the construc- 
tion of a pair (W, W) of exchangeable random variables. Our pair (which 
is somewhat different from those of Stein |St2| and Meckes |Me2| ) satis- 
fies the linearity condition that E(VF / |VK) is proportional to W, and we 
find representation theoretic formulas for quantities such as E(W' — W) k . 
These computations are completely general and apply to arbitrary distri- 
butional approximation. Stein's method of exchangeable pairs is still quite 
undeveloped for continuous distributions other than the normal, but that is 
temporary and as evidence we mention the paper }CFj which develops error 
terms for exponential approximation using quantities like E(W - W) k with 
k small. 

We remark that the bounds in our paper are all given in the Kolmogorov 
metric. Similar results can be proved in the slightly stronger total variation 
metric (see the remarks after Theorem 12. It) . However we prefer to work in 
the Kolmogorov metric as it underscores the similarity with discrete set- 
tings such as |Fuj . where total variation convergence does not occur. We 
also mention that all bounds obtained in this paper are given with explicit 
constants. 



STEIN'S METHOD AND CHARACTERS OF COMPACT LIE GROUPS 



3 



The organization of this paper is as follows. Section [2] gives background 
on Stein's method and normal approximation. Section [3] develops general 
theory for the case that G is a compact Lie group and x T an irreducible 
character. It treats the trace of random elements of 0(n,M), USp(2n,C), 
and U(n,C) as examples. Section [4] extends the methods of Section [3] to 
study spherical functions of compact symmetric spaces. The symmetric 
space setting is natural from the viewpoint of random matrix theory [DnJ , 
[KaS]. After illustrating the technique on the sphere, we treat Dyson's 
circular ensembles as examples, obtaining an error term. 



2. Stein's Method for Normal Approximation 

In this section we briefly review Stein's method for normal approximation, 
using the method of exchangeable pairs [Stlj . For more details, one can 
consult the survey [RR] and the references therein. 

Two random variables W, W are called an exchangeable pair if (W, W) 
has the same distribution as (W',W). As is typical in probability theory, 
let E(^4.|B) denote the expected value of A given B. The following result 
of Stein uses an exchangeable pair (W, W) to prove a central limit theorem 
for W. 

Theorem 2.1. f[Stl]j Let (W, W') be an exchangeable pair of real random 
variables such that E(W 2 ) = 1 and E(W'\W) = (1 - a)W with < a < 1. 
Then for all real xq, 

F(W < x Q ) = / e'^rdx 

V27T J -oo 

< ^gBB +M -L l l -nw>-w\*. 

a V a 

Remarks: 

(1) There are variations of Theorem 12.11 (for instance Theorem 6 of 
[Mel] ) which can be combined with our calculations to prove normal 
approximation in the total variation metric. However Theorem 12.11 
is quite convenient for our purposes. 

(2) In recent work, Roellin [Rlj has given a version of Theorem 12.11 in 
which the exchangeability condition can be replaced by the slightly 
weaker condition that W and W' have the same law. Since ex- 
changeability holds in our examples and may be useful for other 
applications involving Stein's method, we adhere to using Theorem 
I2T1 

To apply Theorem 12.11 one needs bounds on Var(E[(W - W) 2 \W]) and 
K\W — W'\ s . The following lemmas are helpful for this purpose. 

Lemma 2.2. Let (W,W) be an exchangeable pair of random variables such 
that E(W'\W) = (1 - a)W and E(W 2 ) = 1. Then E(W - W) 2 = 2a. 
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Proof. Since W and W have the same distribution, 
E(W'-W) 2 = E(E(W' - Wf\W) 

= E{(W') 2 ) + E(W 2 ) -2E(WE(W'\W)) 
= 2E(W 2 ) - 2E(WE(W'\W)) 
= 2E(W 2 ) - 2(1 - a)E(W 2 ) 
= 2a. 

□ 

Lemma 12.31 is a well known inequality (already used in the monograph 
[Stl| ) and useful because often the right hand side is easier to compute or 
bound than the left hand side. To make this paper as self-contained as 
possible, we include a proof. Here x is an element of the state space X. 

Lemma 2.3. 

Var(E[(W' - W) 2 \W}) < Var{E[(W' - W) 2 \x]). 

Proof. Jensen's inequality states that if g is a convex function, and Z a 
random variable, then g(E(Z)) < E{g{Z)). There is also a conditional ver- 
sion of Jensen's inequality (Section 4.1 of [Du]) which states that for any a 
subalgebra F of the cr-algebra of all subsets of X, 

E(g(E(Z\F))) < E{g(Z)). 

The lemma follows by setting g{t) = t 2 , Z = E((W — W) 2 \x), and letting 
F be the cr-algebra generated by the level sets of W. □ 

3. Compact Lie groups 

This section uses Stein's method to study the distribution of a fixed ir- 
reducible character \ T °f a compact Lie group G. Subsection 13.11 develops 
general theory for the case that x T 1S rea l valued. This is applied to study 
the trace of a random element of USp(2n, C) in Subsection 13.21 and the trace 
of a random orthogonal matrix in Subsections 13.31 and 13.41 Subsection 13.51 
indicates the relevant amendments for the complex setting and Subsection 
13.61 illustrates the theory for U(n,C). 

3.1. General theory (real case). Let G be a compact Lie group and x T 
a non-trivial real-valued irreducible character of G. The random variable of 
interest to us is W = x r (#)> where g is chosen from the Haar measure of 
G. It follows from the orthogonality relations for irreducible characters of 
G that E(W) = and E(W 2 ) = 1. 

The following functional equation will be useful. 

Lemma 3.1. ( / [He2] ; p. 392) Let G be a compact Lie group and \^ an 
irreducible character of G. Then 
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for all a,g € G. 

We now define a pair (W, W) by letting W = x T {g) where g is chosen from 
Haar measure and W = W(ag) where a is chosen uniformly at random from 
a fixed self-inverse conjugacy class of G. Exchangeability of (W, W) follows 
since the conjugacy class of a is self-inverse. Moreover, since x^(o! _1 ) = 
X^(a), one has that x^( a ) 1S rea l f° r an irreducible representations <fi, a fact 
which will be used freely throughout this subsection. 

The remaining results in this subsection show that the exchangeable pair 
(W, W) has desirable properties. 

Lemma 3.2. 

E(W'\W) = ( X ^ a \ ) W. 

\aim(T) J 

Proof. Applying Lemma 13 . 1 1 with <f> = r, one has that 

nW\g) = [ X T (hah- 1 g)dh = ( x ^ g ). 



heG 



\dim(r) 



The result follows since this depends on g only through W. □ 
Lemma 3.3. 

■^-t-K 1 -^)- 

Proof. This is immediate from Lemmas 12.21 and 13.21 □ 

For the remainder of this subsection, if <f> is an irreducible representation 
of G, we let m ( k(T r ) denote the multiplicity of 4> in the r-fold tensor product 
of t (which has character {x T ) r )- 

Lemma 3.4. 

E[(^Tl 9 ]=£^)J^M, 

where the sum is over all irreducible representations of G. 
Proof. Write (W) 2 = Y^<j) m (f>( T ' 2 )x^(9')- Lemma I3TT1 gives that 

nxH(/)\9] = j G XHhah- 1 g)dh = jj^L^fo), 

and the result follows. □ 
Lemma 1331 writes Var([E(W — W) 2 ^]) as a sum of positive quantities. 
Lemma 3.5. 

\ aim{<p) aim(T) J 

where the star signifies that the sum is over all nontrivial irreducible repre- 
sentations of G. 
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Proof. By Lemmas 13.21 and 13. 4| 

E((W'-W) 2 \g) = E[(W') 2 \g] — 2WE(W'\g) + W 



\ dim(<p) dim{r) J 
The orthogonality relation for irreducible characters of G gives that 

Finally, note that 

Var{[E(W' - W) 2 \g\) = E[E({W' - W) 2 \g) 2 ] - (E(W' - W) 2 ) 2 , 

and since the multiplicity of the trivial representation in r 2 is 1, the result 
follows from Lemma 13.31 □ 



Lemma 3.6. Let k be a positive integer. 

(1) e(W - w) k = £i (-i) fc - r C) EWT r W(T*-)i2g). 

(2) E(W - W) A = m,(r 2 ) 2 [8 (l - - 6 (l - 
Proof. For the first assertion, note that 

E [(w' - w) k \g] = ( h ) x T (g) k - r mwr\g}. 

Arguing as in Lemma 13.41 gives that this is equal to 

Thus E(W - W) k is equal to 

E(E[(W' - W) k \g\) 



r=0 

For the second assertion, note by the first assertion that 
E(W - Wf = J>ir (f) £ m^Om^r 4 "') 

r=0 <A 



dim(4>) 



STEIN'S METHOD AND CHARACTERS OF COMPACT LIE GROUPS 

If a is the identity element of G, then W' = W which implies that 



r=0 ^ ' d> 



Thus for general a, 



r=0 

Observe that the r = 0, 4 terms in this sum vanish, since the only contribu- 
tion could come from the trivial representation, which contributes 0. The 
r = 2 term is 

m (r 2 ) 2 . 



^ \ dim{4>) ) 



The r = 1, 3 terms are equal and together contribute 



1 - 



1 - 



dim{r) 



X T (5) 4 



9 GG 



This completes the proof. 

Putting the pieces together, one obtains the following theorem. 



□ 



Theorem 3.7. Let G be a compact Lie group and let t be a non-trivial 
irreducible representation of G whose character is real valued. Fix a non- 
identity element a with the property that a and a -1 are conjugate. Let 
W = X t {q) where g is chosen from the Haar measure of G. Then for all 
real Xo, 

¥(W <x )~ 



x " 

e 2 dx 



< 



\ 



■2\2 



dim(4>) J 



+ 



Here a = 1 — J^tjj the first sum is over all non-trivial irreducible repre- 
sentations of G, and the second sum is over all irreducible representations 
ofG. 



1/4 
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Proof. One applies Theorem 12.11 to the exchangeable pair (W, W) of this 
subsection. By Lemmas 12.31 and 13.51 the first term in Theorem 12.11 gives the 
first term in the theorem. To upper bound the second term in Theorem 12. 1[ 
note by the Cauchy- Schwartz inequality that 



e\w' - w\ 3 < y/E(w - w) 2 e(w - wy. 

Now use Lemma 13.31 and part 2 of Lemma 13.61 □ 

3.2. Example: USp(2n,C). This subsection studies the distribution of 
X T (<?)) where r is the 2n dimensional defining representation of USp(2n, C). 
The only representation theoretic fact needed is Lemma 13.81 which is the 
k = 2 case of a formula from page 200 of [Suj giving the decomposition of r k 
into irreducible representations. In its statement, we let xi, xj , ■ ■ ■ , x n ,x~ l 
denote the eigenvalues of an element of USp(2n, C). 

Lemma 3.8. For n > 2, the square of the defining representation of the 
group USp(2n,C) decomposes in a multiplicity free way as the sum of the 
following three irreducible representations: 

• The trivial representation, with character 1 

• The representation with character x i + x ^ 1 ) 2 + \ + x i~ 2 ) 

• The representation with character Xi+x~ 1 ) 2 — | ^i{x 2 +x~ 2 ) — 
1 

Remark: Lemma [3. 81 could also be easily guessed (and proved) by looking 
at the character formulas for USp(2n,C) on page 219 of [W]. 

Theorem 3.9. Let g be chosen from the Haar measure ofUSp(2n, C), where 
n>2. Let W(g) be the trace of g. Then for all real xq, 



1 f x ° x 2 
F(W < x ) - —= \ e'^dx 
V2it J-00 



n 



to 



1/4 

. Since these bounds hold for all 9 and are continuous in 



Proof. One applies Theorem 13.71 with r the defining representation, and a 
an element of type {xf 1 , • • • , x^ n } where x\ = ■ ■ ■ = x n -\ = 1 and x n = e %e '. 
Then a is conjugate to a~ l and one computes that a = L_££iM_ Using 
Lemma 13.81 one calculates that the first error term in Theorem 13.71 is equal 
to 2 V 4cog (ffL ~^ cos W + ^ _ Q ne computes that the second error term is equal 

24(l-eos(6>)) ~ 
7r(2n+l) 

6, the bounds hold in the limit that 9^0. This gives an upper bound of 
< & as claimed. □ 

3.3. Example: SO(2n + 1, R). We investigate the distribution of x T (9)i 
where r is the 2n + 1-dimensional defining representation of SO(2n + 1,R). 
The only ingredient from representation theory needed is Lemma [3.10l which 
is the k = 2 case of a formula from page 204 of [Su] giving the decomposition 



STEIN'S METHOD AND CHARACTERS OF COMPACT LIE GROUPS 



9 



of r k into irreducible representations (it is also easily obtained by inspect- 
ing the character formulas on page 228 of |Wj). In its statement, we let 
xijX^ 1 , - ■ ■ jin,!" 1 , 1 be the eigenvalues of an element of SO(2n + 1,R). 

Lemma 3.10. For n > 2, the square of the defining representation of 
SO(2n + 1,R) decomposes in a multiplicity free way as the sum of the fol- 
lowing three irreducible representations: 

• The trivial representation, with character 1 

• The representation with character |Q^i Xi+x' 1 )" 2 ' + \ Yli( x l + x i~ 2 ) + 

• The representation with character §QZi x i+ x i 1 ) 2 — \ Yli( x l + x « 2 ) + 

This leads to the following theorem. 

Theorem 3.11. Let g be chosen from the Haar measure of SO(2n + 1,R), 
where n > 2. Let W(g) be the trace of g. Then for all real xq, 



1 f x ° x 2 
¥(W < xo) - \ e~ — dx 

y2vr J-oo 



n 



Proof. One applies Theorem 13.71 taking r to be the defining representation, 
and a to be an element of type {xf 1 , • • • , x^ n , 1} where x\ = ■ ■ ■ = x n -\ = 1 
and x n = e td (i.e. a is a rotation by 6). Then a is conjugate to a~ l and 
a = 2 ^2n+i^ • Using Lemma 13.101 one computes that in the 6 — > limit 
the first error term in Theorem 13.71 is equal to — . One calculates that 

11 | 1/4 

. The proof of the 



the second error term is equal to 



12(2n+l)(l-cos(fl)) 
7rn(2n+3) 



theorem is completed by noting that this goes to as 9 — » 0. □ 



3.4. Example: 0(2n,M). We consider the distribution of x T {g)i where r is 
the 2n-dimensional defining representation of 0(2n,R). The only represen- 
tation theoretic information needed is Lemma 13.121 which is the k = 2 case 
of a result of Proctor [Prj (and also not difficult to obtain from the character 
formulas on page 228 of [Wj). In its statement, we let xijX^ 1 , ■ ■ ■ ,x n ,x~ l 
be the eigenvalues of an element of 0(2n, R). 

Lemma 3.12. For n > 2, the square of the defining representation of 
0(2n, R) decomposes in a multiplicity free way as the sum of the follow- 
ing three irreducible representations: 

• The trivial representation, with character 1 

• The representation with character | (^ Xi + xj) 2 — 4 J2i { x 1 + ^I 2 ) 

• The representation with character \ (^ Xi + ~Xi) 2 +\ J2i { x i + ^T 2 ) — 
1 

This leads to the following result. 
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Theorem 3.13. Let g be chosen from the Haar measure of 0(2n,M), where 
n > 2. Let W(g) be the trace of g. Then for all real xq, 



1 f x ° x 2 
\W < x ) - —j= I e' — dx 

V27T 



CO 



n 



Proof. Apply Theorem 13.71 with r the defining representation of 0(2n, R). 
We take a to be an element of type {xf 1 ,--- , x^ n } where x\ = ■■■ = 
x n -\ = 1 and x n = e %e (i.e. a is a rotation by 9). Then a is conjugate 
to aT l and a = 1 cos ^ , Lemma 13.121 gives the decomposition of r 2 into 
irreducibles, and from this one calculates the first error term in Theorem 
13.71 and sees that in the 9 — » limit it is equal to • ^ ne computes 



that the second error term is equal to 



24n(l-cos(6>)) 
7r(n+l)(2n-l) 



1/4 

. The proof of the 



theorem is completed by noting that this goes to as 9 — » 0. □ 

3.5. General theory (complex case). Let G be a compact Lie group and 
r be an irreducible representation of G whose character is not real valued. 



The random variable of interest to us is W = yx T {g) + X r (<?)J> where 
g is chosen from the Haar measure of G. It follows from the orthogonality 
relations for irreducible characters of G that E(W) = and E(VF 2 ) = 1. 

We now define a pair (W, W) by letting W be as above and W = W(ag) 
where a is chosen uniformly at random from a fixed self-inverse conjugacy 
class of G. As in Subsection 13.11 the pair (W, W') is exchangeable and all 
X^(a) are real. 

The remaining results in this subsection are proved by minor modifications 
of the arguments in Subsection 13.11 



Lemma 3.14. 



Lemma 3.15. 



Lemma 3.16. 



E{W'\W) = ( W. 

\dim{T) J 

E(W' -Wf = l(\- -f^r) 

\ dim{T) J 



^) 2 i 9l = iE™*[(-+^]J^)X*( 9 ). 

where the sum is over all irreducible representations of G. 
Lemma 3.17. 

Var m( W - Wf lg]) ^f mM r + rff (l + - « ) ' , 

where the star signifies that the sum is over all nontrivial irreducible repre- 
sentations of G. 
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Lemma 3.18. Let k be a positive integer. 
(1) E(W - W) k is equal to 

k 



fe/2 



2 k/2 



-1 



,k—r 



r=0 



"V[(r + r ) r W[( r + r 



— \k— ri 



(2) E(W-W) 4 = £W(t + t) 2 ] 2 



2 1 



dim(r) 



x («) 

dim((j)) 
1 



dim(<f>) J 



Finally, one obtains the following central limit theorem. 

Theorem 3.19. Lei G be a compact Lie group and r an irreducible repre- 
sentation of G whose character is not real valued. Let a / 1 fc such that 

a and a -1 



are conjugate. Let W 



i 



X T (g) + X T (g) ) where g is chosen 



from the Haar measure of G. Then for all real xq, 



P(W < x ) 



1 



e 2 dx 

-oo 



< 



5>,[(T + r)*]2 



x'V) \ 

dim(4>) J 







+ 








x r («) 


dim(r) 



3 

2a 



dim(4>) J 



1/4 



Lfere a = 1 — J^t~n , ,/jrs£ siim is ouer a// non-trivial irreducible repre- 
sentations of G, and the second sum is over all irreducible representations 
ofG. 

3.6. Example: U(n, C). Every element of C/(n, C) is conjugate to a diago- 
nal matrix with entries (xi, • • • , x n ) and the representation theory of U(n, C) 
is well understood (see for instance [BuJ). The irreducible representations 
of U(n, C) are parameterized by integer sequences Ai > A2 > • • • > A n . 
The corresponding character value on an element of type {x\,--- ,x n } is 
given by the Schur function s\(x\, • • • , x n ). (The usual definition of Schur 
functions requires that A n > 0, so if A n = —k < 0, this should be inter- 



preted as (a?i 



Cs A+(fc) n ) where A + (k) n is given by adding k to each 



of Ai, • • ■ , A n ). The complex conjugate of a character with data Ai > A2 > 
• • • > A n has data — A n > — A n _i > • • • > — Ai. 

Combining the above information with Theorem 13.191 one obtains the 
following result. 

Theorem 3.20. Let g be chosen from the Haar measure of U{n,C), where 
n>2. Let W(g) = -^[Tr(g) +Tr(g)] where Tr denotes trace. Then for all 
real xq, 



HW < X ) 



e 2 dx 

-00 



< 



n 
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Proof. One applies Theorem 13.191 with r the ra-dimensional defining repre- 
sentation. We take a to be an element of type {x\, • • • , x n } with x\ = ■ ■ ■ = 
Xn-2 = 1, x n -\ = e ld , and x n = e~ ld . Then a and a~ l are conjugate and 
a = 2 ^~T ■ By the Pieri rule for multiplying Schur functions (page 73 
of [Mac] ) . the decomposition of the character of (r + r) 2 in terms of Schur 
functions is given by 



= S (1)S(1) +2— — 



s (2) + + 2 



s (l n ) + s (2,l™- 2 ) 



+ S(-2) + s (-l,-l) 



•2-1 ^-n 

= 2 + S( 2 ) + S (lj i) + S(_ 2 ) + S(-l-l) + 2S( lj0 ™-2 _!). 

One then computes that in the 6^0 limit, the first error term in Theorem 
13. 191 is equal to 2 ^"^}~ 2 < O ne computes that the second error term is 

1/4 

. The proof is completed by noting that this 



equal to 



12(2n-l)(l-cos(6>)) 
vr(n 2 -l) 



approaches as 6 — > 0. □ 

4. Compact Symmetric Spaces 

This section extends the methods of Section [3] to study the distribution of 
a fixed spherical function uj t on a random element of a compact symmetric 
space G/K. Subsection 14.11 gives general theory for the case that uj t is real 
valued. This is illustrated for the sphere in Subsection l4.2l giving a different 
perspective on a result of |DFj and |Melj . We note that since compact Lie 
groups can be viewed as symmetric spaces (see Section 14. 1 f) , the examples 
in Subsections 13.21 13.31 and 13.41 give further examples. Our theorems should 
also prove useful for Jacobi-type ensembles arising from other root systems 
(see for instance f\'r ). 

Subsection 14.31 indicates the changes needed to treat spherical functions 
uj t which are not real valued, and Subsections 14.41 and 14.51 study the trace of 
elements from Dyson's circular orthogonal and circular symplectic ensembles 
as special cases (the circular unitary ensemble is equivalent to U(n, C), so 
was already treated in Subsection I3.6|) . Central limit theorems are known 
for the trace of an element from the circular ensembles (see pages 48-9 of 
[Ra] and [BF]), but our approach gives an error term. 

4.1. General theory (real case). To begin we recall some concepts about 
spherical functions of symmetric spaces. Standard references which contain 
more details are |Hel] . [He2| . [Te| . and [Vj. Chapter 7 of [Mac] is also very 
helpful. 

A Riemannian manifold X is said to be a symmetric space if the geodesic 
symmetry a : X X with center at any point xq is an isometry. Then X 
can be identified with G/K, where G is a connected transitive Lie group of 
isometries of X, and K is a compact group which up to finite index is given 
by K = {g e G : gx = x }. 
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A function E L 2 (G/K) is called spherical if 0^(1) = 1 and the following 
functional equation is satisfied: 

L0 4 ,(xky)dk = u^xju^y) Mx,y G G. 

This equation implies that uj^ is K bi-invariant (i.e. (^^(kigk^) = uj^(g) for 
all ki, k,2 in K), which justifies our writing uJ^g) instead of uj^(gK). 

One reason that spherical functions are important is that if G/K is com- 
pact and Hfi is the G-invariant subspace of L 2 (G/K) generated by u^, then 
is a finite dimensional irreducible representation of G and L 2 (G/K) is a 
direct sum of all such H^. We let dim(4>) denote the dimension of H^. 

In reading this section it is useful to keep in mind that a compact Lie 
group U can be viewed as a compact symmetric space. Indeed, one can take 
G = U x U and K the diagonal subgroup of U; then G/K is identified with 
U under the mapping (u\,U2)K ^> uiu^ 1 ■ The spherical functions uu of 
G/K are indexed by irreducible representations 4> of U and are precisely the 

character ratios moreover, dim{H^) = x^l) 2 - 

Let u T be a non-trivial real valued spherical function of G/K. We are 
interested in the distribution of oo T {g) (normalized to have variance 1). Here 
gK is chosen from the "Haar measure" \i on G/K, i.e. the unique G- 
invariant measure on G/K which satisfies 

/ f(g)Haar{dg) = [ ( [ f(gk)Haar(dk)) n(dgK). 
Jg Jg/k \Jk J 

The following orthogonality relation will be used; see for instance page 45 
of |Klj for a proof. 

Lemma 4.1. 

In particular, Lemma 14.11 implies that W := [dim{r)\ l / 2 uj T has mean 
and variance 1. 

The following lemma is immediate from the functional equation for uu 
and K bi-invariance of u^. 

Lemma 4.2. Let G/K be a compact symmetric space, and a spherical 
function of G / 'K . Then 

I u 4> {kiak 2 g)dkidk 2 = u^aju^g) 

JKxK 

for all a,g G G. 

We define the pair (W,W) by letting W = [dim{T)] l l 2 uj T (g) where gK 
is from the "Haar measure" of G/K, and W 1 = W(ag) where a is chosen 
uniformly at random from a fixed double coset KaK ^ K which satisfies 
the property that KaK = Ka~ l K. Since Ka~ 1 K = (KaK)^ 1 , it follows 
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that (W, W') is exchangeable. Moreover the integral formula for spherical 
functions (page 417 of |He2j ) implies that all ui^a) are real. 

The analysis of the exchangeable pair (W, W) can be carried out exactly 
as in Subsection 13 .14 using Lemmas 14.11 and 14.21 instead of the orthogonality 
relations for compact Lie groups and Lemma 13.11 Hence we simply record 
the results. 

Lemma 4.3. 

E(W'\W) = uj T {a)W. 

Lemma 4.4. 

E(W' -W) 2 = 2(1- io T (a)). 

In the statements of the remaining results, we define the "multiplicity" 
m ( p(T r ) by the expansion 

1/2 

The numbers m^(r r ) are real and non-negative (argue as on page 396 of 
[MacJ with sums replaced by integrals) but need not be integers. Note that 
by Lemma l4.lt 

m4,(T r ) = [dim((j))dim(T) r ] 1/2 / u r {g) r uj^(g). 

Jg/k 

Lemma 4.5. 

E[(W') 2 \gK] = Y j m 4> {T) 2 [dim{<l 1 )} l l 2 uj tj) {a)oj (j> {g), 

where the sum is over all spherical functions ofG/K. 
Lemma 4.6. 

* 

Var([E(W - W) 2 \gK}) = £ m^r 2 ) 2 (1 + u^a) - 2w T (a)) 2 , 

<t> 

where the star signifies that the sum is over all nontrivial spherical functions 
ofG/K. 

Lemma 4.7. Let k be a positive integer. 

(1) E(W> - W) k = Er=o(-l) k ~ r E^K)m (r^H(a). 

(2) E(W> - WY = E m^r 2 ) 2 [8(1 - u T {a)) - 6(1 - ^(a))]. 

Finally, one has the following central limit theorem. 

Theorem 4.8. Let G/K be a compact symmetric space and let uj t be a 
non-trivial real-valued spherical function of G/K. Fix an element a K 



dim((p) 
dim(r) r 
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such that KaK = Ka 1 K. Let W = [dim{r)\ l l 2 ijj T (g) where gK is chosen 
from the "Haar measure" of G/K. Then for all real xq, 

1 r x o 



F(W < Xq) 



e 2 dx 



< 



\ 



+ 



■2\2 



(l-uty(a)) 



-E m ^ 2 ) 2 ( 8 - - c 1 - 

7T V « / 



1/4 



Here a = 1 — w T (a), £/ie /irsi sjim is over all non-trivial spherical functions 
of G/K, and the second sum is over all spherical functions of G/K. 

4.2. Example: the sphere. This subsection studies the unit sphere in 
R n , viewed as the symmetric space SO(n,R)/ SO(n — 1,R). Chapter 9 of 
[V] is a good reference for the spherical functions of this symmetric space, 
and Chapter 4 of |DyM| is a very clear textbook treatment for the special 
case n = 3. Letting ej, • • • ,e n be the standard basis of R n and embedding 
SO(n—l,M) inside SO(n, R) as the subgroup fixing e n , then Kg\K = Kg2K 
if and only if gi{e n ) and 52 have the same last coordinate. From this 
it is not difficult to check that KgK = Kg _1 K for all g, and that the 
double cosets of SO(n — 1,R) in SO(n,M) are parameterized by x n , the 
final coordinate of a point (x\, ■ ■ ■ ,x n ) on the sphere. In what follows we 
let x denote x n . 

From page 461 of [Vj, the spherical functions 101 are parameterized by 
integers I > and satisfy 

Z!(n- 3)1 
y 1 (l + n-3)\ 1 y ' 

Here Cf is the Gegenbauer polynomial, defined by the generating function 

CP{x)t l = {l-2xt + t 2 )-P. 



i>0 



For instance, 
and 



1, C{ (x) = 2px, C 2 p (x) = -p + 2p(l + p)x 1 



LUq(x) = 1, UJ\(x) = X, U2(x) 



nx 



1 



n — 1 

By page 462 of [V], dim(0 = {2l+ X-2W~^ ■ 

We study the random variable W(x) = [dim(l)] 1 / 2 u;i = y/nx. In fact 
sharp (up to constants) normal approximations for W are known: see Dia- 
conis and Freedman [DFj and also Meckes Mel] , who uses Stein's method 
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to obtain an error term of in total variation distance. Our viewpoint 
leads to the following result. 



Theorem 4.9. Let W = y/nx, where x is the last coordinate of a random 
point on the unit sphere in W l . Then for all real xq, 



V(W < xq) 



1 



2 dx 



< 



2V2 



n 



1 



Proof. We apply Theorem 14.81 with r = uj\. Writing uj\{x) 2 as a linear 
combination of uq(x) ando^Oc), one computes that rn^(r 2 ) = 1,TO( 2 )(t 2 ) = 

2> " L ' and that all other multiplicities in r 2 vanish. Letting a be less 



n+2 



than one but close to 1, one computes that the first error term in Theorem 
4.81 is exactly j|^f . Letting a tend to 1 from below, the second error term 
in Theorem H31 vanishes and the result follows. □ 



4.3. General theory (complex case). In this subsection G/K is a com- 
pact symmetric space and tu T is a spherical function which is not real val- 



1/2 



ued. The random variable we study is W = Uil "£' > ^u T (g) + u> T (g) 

As in Subsection 14.11 let W' = W(ag) where a is chosen uniformly at ran- 
dom from a fixed double coset KaK 7^ K which satisfies the property that 
KaK = Ka~ 1 K. Then (W, W) is exchangeable and all uiAa) are real. 

The remaining results in this subsection extend those of Subsection 13. 5| 
and are proved by nearly identical arguments, using Lemmas 14.11 and 14.21 

Lemma 4.10. E(W f \W) = u T (a)W . 

Lemma 4.11. E(W - W) 2 = 2(1 - u T [a)). 

For the remaining results in this subsection, we define the "multiplicity" 
m 4>[( T + ^) r ] by the expansion 

1 1/2 



(uj t + io T y = y 



dim((f>) 
dim{r) 



Arguing as on page 396 of [Mac], one has that the numbers m^[(r + r) r ] 
are real and non-negative (though not necessarily integers). Note that by 
Lemma 14.11 

mj(r + r) r ] = [dim(<j))dim(T) r } 1/2 / (uj T (g) + uj T (g)) ujJg). 

Jg/k v ' 

Lemma 4.12. 

E[(W') 2 \9K] = ~X>*[(r + r) 2 ]^m(0) 1 /2^( a)u;(/) ( 5 ) ) 



where the sum is over all spherical functions of G/K. 
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Lemma 4.13. 

Var([E(W - W) 2 \gK}) = -J^ m ^ T + + ~ ^M)f 



where the star signifies that the sum is over all nontrivial spherical functions 
ofG/K. 

Lemma 4.14. Let k be a positive integer. 
(1) E(W - W) k is equal to 

k 



, — n \ / a. 



2 k/2 

r=0 

(2) K(W - W) 4 is equal to 

— n2t2 



uty(a). 



^m [(r + r) 



2(l-w T (a))--(l-^(a)) 



Putting the pieces together, one has the following central limit theorem. 

Theorem 4.15. Let G/K be a compact symmetric space and let lo t be a 

spherical function of G/K which is not real valued. Fix an element a K 

such that KaK = Ka~ 1 K. Let W = \J ' dim ^ (uj T {g) + uj T (g) S j , where gK 
is from the "Haar measure" of G/K . Then for all real xq, 

1 



P(W < x ) 



x ° 

e 2 dx 

-oo 



< - 



^m [(r + r) 2 p 



1 



- (l-aty(a)) 



+ 



TT * — ' 



T + T 



\2i2 



2a 



(1 -W0(a)) 



1/4 



Here a = 1 — uj T (a), the first sum is over all non-trivial spherical functions 
of G/K, and the second sum is over all spherical functions of G/K. 

4.4. Example: U(n,C)/0(n,~R). This symmetric space can be identified 
with the set of symmetric unitary matrices by the map g i— > gg T (see [Dnj 
or |Fo] for details), and the resulting matrix ensemble is known as Dyson's 
circular orthogonal ensemble. For a thorough discussion of this ensemble, see 
the texts [Foj or [Mtj . In particular, it is known that if a function / depends 
only on the eigenvalues x\, ■ ■ ■ , x n of a matrix from Dyson's ensemble, then 



' G/K 



f=WKI /on,-,**) n 



l<i<j<n 



k=l 
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In this integral, T n is the n-dimensional torus with coordinates 

Si, ; ^ni Xi ^ ^5 \Xi\ — 1' 

It is convenient to let the inner product (/, g) of two functions be defined 
by 

[r(3/2)]» f T ™dx k 

2 > l<i<j<n fc=l 

The spherical functions for this symmetric space are parameterized by 
integer sequences Ai > A2 • • • > A n and are tv\ := p^jjzr. 'i"^ ; the normal- 
ized Jack polynomials with parameter 2. An excellent reference for Jack 
polynomials is Section 6.10 of [Mac] . There one assumes that A n > 0, so 
if A n = —k < 0, P\ should be interpreted as (x± • • • x n )~ k P\ + ( k yi, where 
A + (k) n is given by adding k to each of Ai, ■ ■ ■ , A n . 

To describe some useful combinatorial properties of Jack polynomials, 
we use the notation that if A is a partition and s is a box of A, then 
l'(s),l(s),a(s),a'(s) are respectively the number of squares in the diagram 
of A to the north of s (in the same column), south of s (in the same column), 
east of s (in the same row), and west of s (in the same row). For example 
the box marked s in the partition below 

□ □ □ □ □ 

□ E □ □ □ 

□ □ 

□ □ 

would have l'(s) = 1, l(s) = 2, a'(s) = 1, and a(s) = 3. 

Letting A be a partition of n, and using this notation, two useful formulas 
are the "principal specialization formula" (page 381 of [Mac] ) 



p a (V-- ,i;2) = n 

seA 



n + 2a'(s) -l'(s) 
2a(s) + l(s) + l 



and the formula 

(n + 1 + 2a' (s) - l'{s)) (n + 2a' (s) - l'(s)) 



dim(X) = Y\ 



(2a(s) + l(s) + 2) (2a(s) + l(s) + 1) ' 

A 

which follows from the formula for (P\,P\) on page 383 of [Mac] and the 
fact (Lemma 14. ip that 



(Pa, Pa) 



Theorem 4.16. Let W = |y 1 + \ yTr(g) + Tr(g)j, where g is random 
from Dyson's circular orthogonal ensemble, Tr denotes trace, and n > 2. 
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Then for all real xq, 



1 



e 2 dx 

-oo 



< 



Proof. Apply Theorem 14,151 to the spherical function r = u)^(g) = — -jf^-- 
To compute m^[(r + r) 2 ] for all cj>, one has to decompose (uiru + ^{T)) 2 
into spherical functions, which is equivalent to decomposing (P(i) +-P(i)) 2 in 
terms of Jack polynomials. From the Pieri rule for Jack polynomials ([MacJ, 
page 340), one calculates that 



p (l) p (l) + 2 " 



P (1) P (1 „-i) 



+ p (i) p (i) 



P(2) + gP(i.) + 2 



2n 

i+I r (l") 



P lM +P 



(2,1—3) 



' ' ' %n 



+ P (2 ) + -P( l2 ) 



An 



n + 



J + P (2) + 3^(12) + 2P (1>0 n-a _ 1} + P {2 ) + -P(12). 



We choose a to be an element of type (!,-•• , 1, e ld , e Then KaK 



Koi l K and a 



2(l-cos(9) 



One computes that in the 8^0 limit, the 
first error term of Theorem 14.151 is ^-w ^i^^i^i^^ < ~. The proof is 



completed by computing that the second error term is 
which goes to as 9 — > 0. 



24(n+l) 2 (l-cos(6>)) 
7rn 2 (n+3) 



1/4 
□ 



4.5. Example: U(2n,C)/USp(2n,C). This symmetric space corresponds 
to Dyson's circular symplectic ensemble ([Dn]); see [Fo or [MtJ for back- 
ground on this ensemble. In particular, it is known that if / depends only 
on the eigenvalues x±, ■ ■ ■ , x n of a matrix from this ensemble, then 



Jg/k (2n)! 7 T n ax 



n 



dx k 



l<i<j'<n 



Jl 11 2vr ' 

fe=l 



where T n is as in the previous example. We let the inner product (/, g) of 
two functions be defined by 



(J, 9) 



2" 



l<i<j'<n 



fe=l 



The spherical functions for this symmetric space are parameterized by 

integer sequences Ai > A2 • • • > A n and are lj\ := -7777— ^frf^ , the normal- 

ized Jack polynomials with parameter 1/2. As mentioned in the previous 
example, Jack polynomials are usually defined assuming that A n > 0, so 
if X n = —k < 0, P\ should be interpreted as {x\ ■ ■ ■ x n ) P\+(k) n i where 
A + (k) n is given by adding k to each of Ai, • • • , A n . 
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Letting A be a partition of n and using the notation of the previous 
example, two useful formulas are the "principal specialization formula" (page 
381 of [Mac] ) 



n 



n + 



of (a) 



l'(s) 



and the formula 



(n+S&l- i'( a f) (2n - 1 + a'(s) - 2l'(s)) 
+ Z(s) + l) (a(s) + 2Z(s) + 1) 



s6A 



which follows from the formula for (P\,P\) on page 383 of [Mac] and the 
fact (Lemma 14. ip that 



dim(X) 



(Px,Px) ' 



Theorem 4.17. Lei W = y 1 — ^ [Tr(g) + Tr(g)j, where g is random 
from Dyson's circular symplectic ensemble and n > 2. Then for all real xq, 



1 



e 2 dx 

2ir J-oo 



< 



n 



Tr(g) 



Proof. Apply Theorem 14.151 to the spherical function r = un\ (g) 
To compute m^[(r + r) 2 ] for all </>, one has to decompose (uJm + ^(Tj) 2 
into spherical functions, which is equivalent to decomposing (P(i) + P(i)) 2 in 
terms of Jack polynomials. Prom the Pieri rule for Jack polynomials ([Mac], 
page 340), one calculates that (Pm + P(i)) 2 is equal to 

p d) p (i)+ 2 { h - + p (i) p d) 



P(2) + -P(i*) + 2 



•^1 ' ' ' %n 

2^T P (l n ) + P (2,l""2) 



' ' ' 



2- 

3 



+ ^(1S)+ J P(2) 



2» 



2ra 



J + P (2) + gP(l2) + 2P (li0 n-2 _ 1} + -P(12) + P (2) . 



Now take a to be an element of type (1, • ■ ■ , 1, e lS , e tS ); then KaK = 
KoT^K and a = 2 ( 1 c° s ( s ) _ Q ne com putes that in the # ^ limit, the 



first error term of Theorem 14.151 is _L / 8 ^ n ,. ^, +5 " 3 ) < 1. The second 

" 6(2n-l)(4n-5)(l-cos(fl)) 1 X / 4 
7rn 2 (2n-3) 



error term is computed to be 
8^0. 



which goes to as 
□ 
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